Levy-Lieb constrained-search formulation as a minimization of 

the correlation functional. 

Luigi Delle Site 

Max-Planck-Institute for Polymer Research 
Ackermannweg 10, D 55021 Mainz Germany. 

Abstract 

The constrained-search formulation of Levy and Lieb, which formally defines the exact 
Hohenberg-Kohn functional for any iV-representable electron density, is here shown to be equivalent 
to the minimization of the correlation functional with respect to the N — 1 conditional probability 
density, where N is number of electrons of the system. The consequences and implications of such 
a result are here analyzed and discussed via a practical example. 
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I. INTRODUCTION 



The Hohenberg-Kohn (HK) theorem [l| has opened new perspectives to the calculations 
of electronic-based properties of condensed matter [2|, and, an aspect often disregarded, has 
given profound new insights into the general understanding of quantum mechanics. In fact 
the 3iV-dimensional Schrodinger problem for the ground state of an electronic system: 

H n iP(t u ...r 2 ) = E ^t x , ...r 2 ); H N = £ (~^f) + £ f (rA + £ — (1) 

i=l,N 1 i=l,N i<j r *i 

where v (r«) is the external potential, Di<j — is the electron-electron Coulomb term, E 



is the energy of the ground state and V( r ij ••• r iv) is the 3A^-dimensional antisymmet- 
ric ground state wavefunction 3], is transformed into a "manageable" variational prob- 
lem in three dimensions where the central role is played by the electron density: p(r) = 

iV J f2jv _ i -^*(r, r 2 , r N )ip(r, r 2 , r N )dr 2 ....dr N , where fijv-i is the AT — 1 spatial domain. 

In explicit terms the variational problems is written as: 

E = Min p E[p] (2) 

where f n p(r)dr — N (ft being the spatial domain of definition) and E[p] = T[p] + V ee [p] + 
V ext [p] is the energy functional composed respectively by the kinetic, electron-electron po- 
tential and the external potential functional. However in its original formulation the HK 
theorem and the related variational problem have got a restricted field of applicability; it 
is valid only if the electron density p(r) is f-representable, that is if p(r) is the density cor- 
responding to an antisymmetric wavefunction of the ground-state of an Hamiltonian of the 
form of Eq{T] It follows that the correct formulation of the variational problem becomes: 

E = Min p E v \p] (3) 

where v refers to the u-representability of p(r). As discussed in Ref. j^J, there are no general 
conditions for a density to be f-representable and this makes the use of the HK theorem and 
its associated variational principle not practical. A generalization of the HK theorem which 

n 

does not require p(r) to be f-representable was found, in parallel, by M.Levy [4] and E.Lieb 
j^l and it is usually known as the Levy constrained-search formulation or Levy-Lieb 
constrained-search formulation (j); in this paper we adopt the latter terminology. We 
also notice that recently P.Ayers [7| has further clarified this concept and developed an 



axiomatic treatment of the Hohenberg-Kohn functional. In the following we briefly describe 
the crucial aspects of the abovementioned approach which are relevant for the current work. 
The starting point of the theory is the distinction between the ground state wavefunction, if}, 
and a wavefunction ip x that also integrates to the ground state electron density p(r) . Since 
if) is the ground state wavefunction, we have: 

(if> x \H N \if> x )>(i(>\H N \i{>) = E . (4) 

Taking into account that V^ x t[p] is a functional of p only, EqHJcan be written as: 

(■0A \T + Veel > ^|T + y ee |V) (5) 

where T and V ee are respectively the kinetic and Coulomb electron-electron operator as 
defined in EqJTJ The meaning of EqJS] is that ip is the wavefunction that minimize the 
kinetic plus the electron-electron repulsion energy and integrates to p. It follows that the 
initial variational problem of Eqj2] can be transformed in a double hierarchical minimization 
procedure which formally allows for searching among all the p's which are iV-representable, 
i.e. it can be obtained from some antisymmetric wavefunction; this is a condition which 
is much weaker and more controllable than the w-representability. In explicit terms such a 
formulation is written as: 



E = Min p 



(i/> x \T + V ee \i> x ) + J v(r)p(r)dr 



(6) 



The inner minimization is restricted to all wave functions ip x leading to p(r), while the outer 
minimization searches over all the p's which integrate to N. The original HK formulation can 
then be seen as a part of this new one once its universal functional, F[p] = (ip\T + V ee \ ip) 
is written as: 

F[p] = Min^p (ip \T + V ee \ip) . (7) 

The purpose of this work is to show that F[p] can be determined solely by a minimiza- 
tion with respect to the N — 1 conditional probability density of the electron correlation 
functional. This latter will be shown to be composed by the non local Fisher information 
functional and the electron-electron two-particle Coulomb term. The advantage of this 
representation is manifold; it further clarifies the connection of electronic properties to the 
Fisher theory and shows that the knowledge of such a functional is the crucial ingredi- 
ent in density functional based approaches; it also identifies the Weizsacker kinetic term, 



as necessary component of the universal functional F[p] and, in practical terms, 
offers an objective criterion of evaluation of "approximate" exchange and correlation func- 
tional, i.e. among two functionals, the physically better founded is the "smaller" one. In 
order to show the practical aspects of our idea we illustrate a potential application. 



II. THE NEW REPRESENTATION 



Before writing the functional in the conditional probability density formalism, we need 
to define such a quantity. Let us consider a generic fermionic wavefunction ^(r 1; ....rjy), for 
simplicity we consider a real wavefunction, but the extension to a complex one can be also 
done we do not consider the spin dependence explicitly, however this will not influence 



the main conclusions. Then the iV-particle probability density is 10 



iVV*(ri, ....rjvMn, ....r N ) = 6 2 (n, r N ) (8) 



and this can formally decomposed as [10, lll[ : 



2 (r 1 ,....,r JV )=p(r 1 )/(r 2 , (9) 

where p(i"i) is the one particle probability density (normalized to N) and /(r2, , rjv/ri) 

is the N — 1 electron conditional (w.r.t. ri) probability density, i.e. the probability density 

of finding an iV-1 electron configuration, C(r 2 , , rjv), for a given fixed value of The 

function / satisfies the following properties: 

(«) / /(r 2 , ,rjv/ri)dr 2 dr N = l\/r 1 

(ii) /(n, ..r i ...r i _i,r i+a ...,r J y/r i ) = 0; for i = j; Vi, j = 1, N 
(Hi) f(ri, rj..,r j ...r fe _ 1 ,r fe+ i..,r iV /r fe ) = 0; for i = j;Vi,j^k (10) 

The property (Hi) of Eq{10] assures us that / reflects the fermionic character of an electronic 
wavefunction. In fact it says that if any two particles are in the same 'state' "r" the 
probability of that specific global configuration is zero. In principle, together with condition 
(ii), this is a way to mimic the antisymmetric character of the fermionic wavefunction since 
for fermions | -0 (r i , ...rj, ...rj, ...r N )\ 2 = 0; for i = j, Vi, j.lt must be noticed that condition 
(in) is complementary to (ii). With this formalism the term (ip\T + V ee \ ip) can be written 
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as (see Refs. 



ivpwi 2 , 1 / 



|V r /(r',....,rjv/r) 



fijv-i /(r, ,rjv/r) 

(N-l) lp(r) 



-dr ....dr 



N 



/(r '' > rN/r) -dr....dr N 



r — r 



dr + 
411) 



where we have identified r x with r and made use of the property of electron indistinguisha- 
bility, thus r could be identified with any of the (and the same for r' identified here with 
r 2 ) without changing the results; a further consequence is that the Coulomb expression (last 
term on the r.h.s.) is written as the sum of iV — 1 identical terms for the generic r and r' 
particles. Using EqJTTJthe Levy-Lieb constrained-search formulation can then be written as: 

|Vp(r)' 2 



E = Min p (linj (r[/,p]) + ~ J dr + J v(r)p(r)dr 



where 



r[/,p] = g/p(r 



|V r /(r, ....,r N /r) 



(12) 



fitf-i f(r, ,t n /t) 

(N-l) I p(r) 



-dr ....dr 



N 



/(r ' ^ N,r) -dr....dr N 



dr + 
dr. 



r — r 



(13) 



In this way we have transferred the problem from from ip to / which means that the focus 
is now on T[f,p], i.e., as discussed in Ref. 12|, the correlation functional. 



III. A PRACTICAL EXAMPLE: THE PARAMETRIC EXPONENTIAL FORM 
OF / 



In our previous work [12 ]. we have proposed an approximation for / based on a two- 
particle factorization: 



where 



~E(r) 



e -E H (r, ri ) d 



(14) 



(15) 



, N is the number of particle, and u the volume corresponding to 



here E H (r,Ti) 

one particle. Such an approximation, due to its simplicity, allows us to write an analytic 
expression of the Fisher functional which can be used in a straightforward way in numerical 



5 



calculations. However it does not satisfy the condition (ra) of EqJTQl and, for this reason, 
in order to use it into the Levy-Lieb constrained-search scheme it must be extended. The 
expression we propose here is the following: 

/(r 2 , ...v N /r) = n n=2 ^e f W-^( r ' r ") x H l>m e^ E «^ (16) 

with: 

e -I« = J U n=2:N U l>m e-^^ E «^dr 2 dr N (17) 

Here 7 and (3 are two free parameters. As it can be easily verified this expression of / 
satisfies all the requirements of EqJTOl The meaning of / as expressed in EqJT6] is that the 
probability of finding a certain configuration for the N — 1 particles, having fixed particle 
ri = r, depends not only on the fixed particle and its interaction with the N — 1 other 
particles as before, but also on the mutual arrangements of the N — 1 particles (it has 
also to be kept in mind that using the particle indistinguishability the formalism can be 
applied to any r, as a fixed particle). The parameters 7 and (3 express how important 
the N — 1 mutual interactions are with respect to the interactions with r. Being now / a 
biparametric function, one can use the Levy-Lieb constrained search in our formulation and 
find the optimal values for 7 and (3. This practical example shows two different aspects of 
our formulation; basically we have shown that indeed it is possible to build a function / and 
actually it can be chosen in a way that its optimal expression can be determined via the 
constrained-search formulation. It must be noticed that this form of / is still rather simple 
since the spins are not explicitly considered when constructing the function and thus one 
cannot distinguish between the exchange and the correlation part of the electron-electron 
interaction as it is done in standard Density Functional Theory; as a consequence one should 
expect only an overall average description of these two terms which are here incorporated into 
the global correlation. However the construction of a more complete expression of /, which 
takes care of the effects of the spins, is the subject of current investigation. This emphasizes 
once more the merit of the general procedure shown here, that is different expressions of /, 
with different degrees of complexity, can be proposed and their relative validity checked by 
the constrained-search procedure. 
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IV. DISCUSSION AND CONCLUSIONS 



As anticipated in the introduction, the consequences of Eqs JT2]TT3l are rather interesting. 
The Levy-Lieb variational principle can be reformulated as: The universal functional F[p] 
is the one with the minimum correlation functional with respect to the electron conditional 
probability density. This new interpretation of the HK universal functional tells us that 
only an accurate description of the correlation effects, considering the Weizsacker term as 
a necessary term, leads to an accurate description of the whole energy functional; such a 
criterion is necessary and sufficient. It is obvious that it is necessary; without knowing 
T[f,p], F[p] cannot be known; it is sufficient because once T[f,p] or better f(r 2 , ...r N /ri) 
is (in principle) known than the whole energy functional is known explicitely. Clearly, the 



"true" /(r2, ...rjv/ri) is very difficult if not impossible to obtain 



1311. however it can be 



sufficiently well described on the basis of mathematical requirements and physical intuition 



as done for example in Ref. 



12| and as shown in the previous section. From this point of 



view, Eqs JT2lfT3l can be seen as an objective criterion to design, on the basis of physical 
intuition and fundamental mathematical requirements, valid energy functionals. In fact, as 
done in Ref.jljl and in the previous section, one can construct well-founded expressions 
for / keeping in mind the physical meaning of the electron correlation effects and the 
necessary related mathematical prescriptions of EqJTOl Next one can make use of Eqs JT2][T3l 
and choose among different functional forms of /, the one giving the "smaller" T. It 
must be noticed that in this work we do not claim that finding a functional form of / is 
easier or more rigorous than to find an exchange-correlation functional in standard Density 
Functional Theory; it represents an alternative or complementary approach to the latter. 
However, the approach based on / allows one to express in a more direct way, via the choice 
of different forms of /, the physical principles related to the electron correlation effects and 
to have an explicit form of the correlation term for the kinetic functional which is of great 
advantage for kinetic functional based methods (see e.g. Refs. fl^. fislj ) An important aspect 
linked to the statement above is that the term: I f p(r) f |Vr/(r ,....,r N /r)\ fa' ....dr jydr, is 
the well known non loca 
is available (see e.g. 
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16 



Fisher information functional about which a vast literature 
171 ] and references therein); this term is very often linked to 
the electron correlation functional and electronic properties (see Refs. jla. fl9]). our work 
further clarifies this connection, suggesting that the results known from the analysis of 
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the Fisher functional could be employed in this context. In conclusion we have shown an 
alternative view of the Levy-Lieb constrained search approach and provided an example 
which clarifies the practical advantage of our idea; in this sense the present work it is not 
merely a marginal new formal contribution to a rather well-known method, but gives a new 
powerful insight into the field of applicability for realistic systems. 
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